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Abstract. We consider a hybrid system consisting of two normal metal leads weakly 
connected to a superconductor. Current-current correlations of the normal leads are studied 
in the tunneling limit at subgap voltages and temperatures. We find that only two processes 
contribute to the cross-correlation: the crossed Andreev reflection (emission of electrons in 
different leads) and the elastic cotunneling. Both processes are possible due to the finite size of 
the Cooper pair. Noise measurements can thus be used to probe directly the superconducting 
wave function without the drawbacks appearing in average current measurements where the 
current is dominated by direct Andreev reflection processes. By tuning the voltages it is possible 
to change the sign of the cross correlation. Quantitative predictions are presented both in the 
diffusive and ballistic regimes. 



Introduction. - The Andreev reflection is the process of elastic transfer of two electrons 
forming a Cooper pair from the superconductor to the normal metal £Q. Electrons that are 
emitted in the normal metal cross the normal/superconductor interface within a distance of the 
order of the size of the Cooper pair. It is thus possible for the two electrons to be transmitted 
in two different normal leads, when the distance between the contacts is comparable with 
the pairs' size. This process, called crossed Andreev reflection (CAR), can be used to probe 
directly the spatial structure of the Cooper pair, but how can it be detected? Different 
authors have proposed to measure the conductance matrix in the multiterminal structure and 
to extract the distance dependent contribution of the CAR 2,3,4,5 . However this procedure 
has two main drawbacks. First, CAR comes along with elastic cotunneling (EC) 0, the 
transfer of an electron from one normal lead to the other one through the superconductor 
(conserving its spin and energy). Second, CAR and EC contributions to the average currents 
are dominated by direct Andreev reflection in each normal lead, i.e., by the current associated 
with the transfer of two electrons in the same lead We will show that one can pick up 
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Fig. 1 - Schematic picture of the three terminal device. 

the CAR or the EC contribution directly by measuring the cross-correlation of currents in the 
two normal leads. 

As a matter of fact, due to the discrete nature of charges, current fluctuates in time around 
its mean value. Current noise measurement revealed a powerful probe for electronic systems 
[7j- In particular, cross-correlation in multitcrminal superconducting hybrid structures has 
deserved some attention since it was first anticipated [SJIHl and then predicted [TniimiT^llT^l 
t° snow a s ig n change with respect to normal metallic structures. A simplified 
explanation of this effect is the following. If electrons in the two leads are emitted from 
one Cooper pair, one expects a positive correlation, since both electrons appear at the same 
time in each lead. By contrast, electrons from a fermionic source due to Pauli principle 
arrive one by one in each channel and they are transmitted either in one or the other lead, 
leading to a negative correlation. Most calculations consider Y-shaped structures attached to 
a superconducting lead through a sing le junction |10I11I12I13I15I16I17| . This geometry is the 
standard one for interference in optical experiments, but it leads to a variety of elementary 
processes for charge transfer in electronic samples. This makes the interpretation of the cross- 
correlation's sign more subtle PI - 

In this paper we consider the different situation (depicted in Fig. ^) where each normal 
arm is directly connected to the superconductor through a tunnel interface. We find that 
CAR contributes to the cross-correlation with a positive sign, while EC contributes with a 
negative sign, and there are no spurious contributions from direct Andreev reflection. The 
CAR and EC contributions can be selected independently by tuning the voltages of the normal 
terminals. Cross-correlation in tunneling systems thus provides a direct way (i) to probe both 
CAR and EC and («) to measure the sign change of the correlations. 

We present quantitative predictions for the dependence of the cross-correlation spectral 
noise on the distance between normal arms both in the ballistic and diffusive regimes. 

Effective Hamiltonian. - Let us consider a conventional Bardeen-Cooper-Schrieffer su- 
perconductor weakly connected to two normal metal leads called A and B (see Fig. To 
describe charge transport we use the standard tunnel Hamiltonian: H = Hs + + i?jy + 
Hjt + H^ where H 1 ^ and Hs describe the clean or impurity-disordered normal lead a (with 
a = A or B) and the superconductor, respectively. Explicitly Hfi = J2ka ^fc c fcl c fc<r an d 
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H s = J2qa(q d q<7 d q<T + H q [^ d qi d - q ] + ^* d -q] d ql] with A the superconducting gap. The 
operators c ka and d qa are the destruction operators for electrons of spin a in the normal 
metal a and in the superconductor, respectively. The indices fe and q together with the spin 
indicate quantum numbers labelling the eigenstates in the disconnected leads. The tunneling 
part of the Hamiltonian is given by: 



Ht — 5Z^fcq C fci^qT + tkq d qa C k(j] J 
kqa 



(1) 



where t kq = J drdr' t a (r,r')ipk(r)ip*(r') and t a (r,r') is the quantum amplitude for an 
electron to tunnel from position r in the normal lead a, to position r' in the superconductor. 
We introduce the electrostatic potential by the standard gauge transformation — ► = 
6c + eV a . 

We are interested in the subgap regime defined by \eV A \, |eVs|, ksT <C A. Then, no 
excitation can be created in the superconductor over times larger than h/A. Under these 
conditions, in lowest order of perturbation theory, only two different processes contribute to 
charge transport: Andreev reflection and elastic cotunneling. In both cases the corresponding 
quantum amplitude can be obtained using second order perturbation theory JU] . Specifically, 
for the Andreev process, one can calculate the quantum amplitude for destroying a Cooper 
pair in the superconductor and create two electrons, one in the lead a with spin f and the 
other in lead (3 with spin |. The result is: 



.4 



a/3 
kp 



tkqtp- q u q v q 



tZ+Eq 



4 



E a 



(2) 



where u q = 1 — v q = (1 + ( q /E q )/2 and E q = y A 2 + Q q is the superconducting spectrum. 

The amplitudes A AA and A BB correspond to the injection of both electrons in the same lead, 
while A AB and A BA give the amplitude of the crossed Andreev reflection. 

In a similar way one can calculate the amplitude for clastic cotunneling (B, a) — > [A, a): 



T — + A + B * 

J-kp — b kq L pq 

q 



Eq+HZ 



E -£ p 

^q ?p 



(3) 



It is convenient to write an effective Hamiltonian that takes into account these two pro- 
cesses 1251: 



H cfi = H A 



H B + J2[J a p + h.c.]+T + T^ 



(4) 



a,/3 



where J a p = Y.k P ^tp c k\ c p\ and T = J2k P a T kp c tl c p<y This new Hamiltonian is equivalent 
to the initial one at low energies and it allows to obtain current and noise straightforwardly. 



Currents and cross- correlation. 



Current operator in each normal lead is defined, as 



usual, by the time derivative of the particle number operator N a — ^2 ka c^ a c' k(y 



e &N a /&t = —{ei/K)\H , N a ] (sign conventions are defined in Fig. 2J. We thus obtain: 



26 %C 

■ [2J AA + J AB + Jba +T} + h.c. ; I B = — [2J BB 



Jab + J B A-T] + h.c. (5) 



n n 
It is now possible to evaluate the current in both normal leads and their zero-frequency 



self- and cross-correlation defined as: 



/+oo 
dt([5I a (t),SI (O)}+), 
-OO 



(6) 
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where SI a (t) = I a (t) — (I a ), ol and /? take the values A or B, and the brackets denote quantum 
and statistical averaging. One can then apply standard time-dependent perturbation theory 
to the lowest non vanishing order (~ t 4 ). We obtain: 

(j A ) = ^[{Nt A -Ni A ) + {Ni B -Ni B ) + {Nt c -Nt C )} , (7) 
Saa = ^ [2 (A^ 4 + A^ A ) + + + ( N ec + n ec^ (8) 

Ap 2 

Sab = ^ [(iV^ B + Nt B ) - (N EC + N* c )] , (9) 

with 

= 2whJ2 \Af p \ 2 H^ d (( k , SjSfa + £ p + eV a + eV ) , (10) 

kp 

N EC = 2irhJ2 \T kp \ 2 HI?(Zk,ZpMk - £ P + eV A - eV B ) , (11) 

kp 

and Ht nd (^,c P ) = /(&)/&») , Ht nd (Zk,z P ) = [i - mm - m P )} , #i ot (a,£ P ) = 

/(&)[! " /Up)] = ^ ot (e p ,a), where /(&) = (c^cgj is the Fermi function. 

Eq. (J7J) for the current agrees with the result obtained by Falci et al. in Ref. As 
mentioned above the average current has three components: the direct Andreev (N AA or 
N BB ), the crossed Andreev (N AB ), and the elastic cotunneling (N EC ) currents. The direct 
Andreev term is instead absent in the cross-correlation given by Eq. ©: cross-correlation is 
thus a direct measure of CAR and EC. It is convenient to introduce S CAR — 4e 2 /H 2 (N AB + 
N AB ) and S EC = Ae 2 /h 2 {N EC + N BC ) in terms of which we have: S AB = S CAR - S EC . 
Since N^> B and are always positive, this expression implies that CAR and EC processes 
contribute to Sab with positive and negative sign, respectively. A simple interpretation of 
this fact is that CAR implies instantaneous currents of the same sign in both leads, while EC 
implies instantaneous currents of opposite signs. 

To obtain more quantitative predictions we follow the procedure of Refs. (191120) . The 
main result is given in Eq. 1 )16(1 below, with the amplitudes given in Eqs. 1)19( 1 and 1)22)) for 
the ballistic and diffusive regimes, respectively. 

Introducing U = (V A + V B )/2 and Gq = 2e 2 /h, S CAR can be rewritten as 

S CAR = 2G Q J de A(e) [f(e + eU) (1 - f(e - eU)) + f(e - eU) (1 - f(e + eU))] , (12) 
with F(C;U,tB) = 2iru(()v(0{(E(() + eV A - U)' 1 + (E(() + e^+^r 1 }, and 

A(e) = J d(dCF((; e-eU,e + eU)F(g; e-eU,e + eU) S(e - eU, e + eU; C, CO • (13) 
The function 5 reads: 

S(£a,£b;CO = / dr x dr 2 f dr 3 d ri f dr[...dr' A t A *(r 1 ,r[)t A (r 2 ,r' 2 )t B *(r 3> r' 3 )t B (r A ,r' 4 ) 
JVa JVb JVs 

K U (r x ,r 2 )K iB (r 3) r^K^, r' 3 )K ( , (r' 4 , r[) (14) 



where K^{r\,r 2 ) — J2k H^. ~ ^k)tpk(fi)tp k (r 2 ) is the single particle spectral function, Va, Vb, 
Vs indicate the volumes of leads A and B, and of the superconductor. 
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We assume that tunneling only occurs with constant amplitude and between neighboring 
points lying on the junctions surface S a . Following Ref. 6, 19] one can then trade the t, K^ A , 
and K^ B dependence for the normal-state tunnel conductances per unit surface, g A and g B , 
(plus the density of states of the superconductor, vs): 

3(a,fr»;C, C')-^(C,C) = |yn / d ^ / d 2 r B K ( (r A ,r B )K c ,(r B ,r A ) . (15) 

16ir 2 e 4 is s J Sa J Sb 

Fermi functions in Eq. I(12|l restrict the range of variation of |e| to be at most of the order of 
max(k B T, \eVj\\, \e.V B \) <C A. One can then safely discard voltage and e dependence of F in 
Eq. ((13(1 . A similar procedure can be applied to S EC as well. We finally obtain: 



Sab = 2eGQ 



eV A + eV B \ aC ab. n , w n - ( eV A ~ eV B > EC 



(V A + Vs) coth ( I A— - (V A - V B ) coth I -^-^ ) A 



(16) 

with A i = J dCdC'S'tO^'fOSCCC), where i = {CAi?,£C}, $ C ^(C) = 2vrA/(A 2 + C 2 ) 
and $ (£) = 27r£/(A 2 + ^ 2 ). Eq. l|ltj|) singles out the voltage and temperature dependence 
of the cross-correlation. The spatial dependence due to the coherent propagation of electrons 
in the superconductor between the two junctions is isolated in the amplitudes A CAR and A EC . 
The same amplitudes appear in the expressions for the current obtained from |7|>: 

{I A ) = I° A (V A ) + G Q (V A + V B )A CAR + G Q (V a - V B )A EC . (17) 

Here I A (V A ) is the direct Andreev contribution discussed by Hekking and Nazarov ^Hj giving 
rise to the non linear I-V characteristic. The second and third term on the rhs of Eq. 1)17(1 
give instead the CAR and EC contribution to the current and they have been discussed in 
Refs. |2|31l3lEj- F° r completeness we give the current auto-correlation (noise): 

S AA = 4e/l(Vk)coth (eV A /k B T) + S CAR + S EC . (18) 

The first term of the rhs was obtained in Ref. [20] and dominates the noise. 

Comparing Eqs. I(16|) and 1(17(1 the advantage of measuring the low-temperature cross 
correlation over the mean current is apparent. Indeed, in Sab for T = there is no direct 
term. This allows to measure directly A CAR or A EC by setting Va = V B or Va = —V B , 
respectively. In contrast, the same procedure would not work by measuring the current. As 
we will see below, A CAR fa A EC = A, thus the current becomes (I A ) = I° A {V A ) + 2G Q AV A , 
with no residual dependence on V B . The interesting A dependence is thus hidden by the direct 
term 1° . In order to raise the amplitude degeneracy some authors [T^| fH|l21|ll41i22ll23j proposed 
to perform the experiment with spin-polarized normal leads, and measure the drag current 
Ia(Va = 0, V B ). This is not necessary if S AB is measured. A last remark on the temperature 
dependence is in order. For temperatures k B T ^> | eV^| , \eV B \ cross-correlation is completely 
suppressed. 

We discuss now briefly the amplitude in the two interesting regimes: ballistic and diffusive. 
In the ballistic regime one finds 0J|S]: 



A EC ) ~ -Gl J SA J Sb d rAd rB { co S Hk F \r A -r B \) ) {k F \r A -r B \)> ' (19) 



where ^ = hvp/A is the superconducting coherence length and vp, kp the Fermi velocity 
and momentum. In multichannel junctions, the trigonometric functions in 119(1 are averaged 
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Fig. 2 - Amplitude Ad for two square contacts on a diffusive superconductive film. From top to 
bottom: L/t; = oo, 5, 2, 1, 0.5. Dependence on the distance R for different side lengths L. Inset: 
dependence on the side length for R = 0. 

out, giving the same numerical factor. Thus A CAR ss A EC = A^. For junctions with a typical 
size much smaller than the distance, R, between the two contacts, the spatial dependence is 
given by At ~ e~ 2R '^ b / (kpR) 2 . The exponential decay is related to the characteristic size 
of the Cooper pairs in the superconducting lead. The algebraic prefactor further reduces the 
intensity of the effect on the scale of the Fermi wavelength as discussed in Refs. |4I5I24| . Cross 
contributions are expected to be small in ballistic systems. 

In diffusive superconductors the coherence length £<j = \JhD j A (D is the diffusion con- 
stant) and the distance between the two junctions largely exceed the elastic mean free path 
l e . Impurity averaging is performed on Ijl5(l 22 ,23 within the usual approximation fcpZ e -C 1: 

(K c (r A ,r B )K c (r B ,r A )) imp ='^[P c _ c (r A ,r B ) + P c _ c (r B ,r A )} . (20) 
The Cooperon P e satisfies the diffusion equation in the superconducting lead: 

-hDA ri P c - ieP e = 5 3 (n - r 2 ) . (21) 

Then, we find 

A d = A CAR = A EC = ^f- [ d 2 r A d 2 r B P 2 , A (\r A -r B \). (22) 

If the superconducting electrode is a thin film with a thickness d <§; £<j (characterized by 
its sheet conductance G\j = 2e 2 vsDd in the normal phase), we can use the picture of two- 
dimensional electron diffusion. Then we find P2i/\{FL) — (l/TrhDd)Ko(V2R/ £,d) , where Kq is 
the modified Bessel function of the second order. In a thicker film, P 2 iA (R) = (1/ ' ^hDR)e^ 2R ^ d . 
We calculate explicitly the value of Ad when the normal leads are square films of side L de- 
posited at distance R over a superconducting film. The results are presented in Fig. As 
expected, the distance dependence is essentially determined by P2ia(-R)- The size dependence 
is more interesting. For large values of Ad(R = 0) grows linearly with L, this is due 

to the fact that the main contribution comes from a small stripe of width ^ and length L 
in each contact. This can be used experimentally to increase the signal, there are no draw 
backs in increasing the lateral side of the contacts. For small contact sizes L <C £d, the whole 
contact contributes with the same value giving A d ~ L 4 P 2 iA(R)- 
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Conclusion. - We calculated the current-current cross-correlation in multiterminal hy- 
brid structures consisting of two normal metallic leads attached to a superconductor through 
tunnel barriers. We found that varying the voltage biases between the normal leads and the 
superconductor allows to probe separately two non-local mechanisms for the charge transfer 
in such systems. We evaluated this effect for realistic devices and we found that it may be 
observable in diffusive structures, provided that the two normal arms are deposited at a dis- 
tance not larger than few tens of nanometers (of the order of the superconducting coherence 
length in the diffusive regime). In view of the recent observation of the doubled shot noise in 
a superconducting/normal metal tunnel junction 125| we think that measuring such effects is 
within the reach of present technology. 

* * * 

We thank D. Feinberg for pointing out the advantage of special geometries of contacts. G.B, 
F.P., and F.W.J.H. acknowledge financial support from CNRS/ATIP-JC 2002 and Institut 
Universitaire de France. MH acknowledges financial support from the ACFJC no 2036 from 
the French Minstry of Research. 

REFERENCES 

[I] Andreev A. F., Sov. Phys. JETP, 19 (1964) 1228 (Andreev A. F., Zh. Eksp. Teor. Fix., 46 
(1964) 1823). 

[2] Byers J. M. and Flatte M. E., Phys. Rev. Lett, 74 (1995) 306. 

[3] Deutscher G. and Feinberg D., Appl. Phys. Lett, 76 (2000) 487. 

[4] Recher P., Sukhoroukov E. V. and Loss D., Phys. Rev. B, 63 (2001) 165314. 

[5] Falci G., Feinberg D. and Hekking F. W. J., Europhys. Lett, 54 (2001) 255. 

[6] Averin D.V. and Nazarov Yu. V., Phys. Rev. Lett, 65 (1990) 2446. 

[7] For a recent review see Blanter Y. M., and Buttiker M., Phys. Rep., 336 (2000) 1. 

[8] Anantram M. P. and Datta S., Phys. Rev. B, 53 (1996) 16 390. 

[9] Martin T., Phys. Lett A, 220 (1996) 137. 

[10] Torres J. and Martin T., Eur. Phys. J. B, 12 (1999) 319. 

[II] Torres J., Martin T. and Lesovik G. B., Phys. Rev. B, 63 (2001) 134517. 
[12] Lesovik G. B., Martin T. and Blatter G., Eur. Phys. J. B, 24 (2001) 287. 
[13] SHECHTER M., Imry J. and Levinson Y., Phys. Rev. B, 24 (2001) 287. 

[14] Taddei F. and Fazio R., Phys. Rev. B, 65 (2002) 134522. 
[15] Borlin J., Belzig W. and Bruder C, Phys. Rev. Lett., 88 (2002) 197001. 
[16] Samuelson P. and Buttiker M., Phys. Rev. Lett, 89 (2002) 046601. 
[17] Samuelson P. and Buttiker M., Phys. Rev. B, 66 (2002) 201306. 

[18] For general reviews see also Quantum noise in mesoscopic systems, edited by Nazarov Yu. V. 

(Kluwer, Dordrecht) 2002. 
[19] Hekking F. W. J. and Nazarov Y. V., Phys. Rev. Lett, 71 (1994) 1625 and Phys. Rev. B, 

49 (1994) 6847. 

[20] Pistolesi F., Bignon G. and Hekking F. W. J., |ArXiv:cond-mat /0303165 
[21] Melin R. and Feinberg D., Eur. Phys. J. B., 26 (2002) 101. 

[22] Chtchelkatchev N. M. and Burmistrov I. S.,|ArXiv:cond- mat/0303014 Chtchelkatchev 

N. M. and Ma r'enko M., |ArXiv:cond-m at/0306552| 
[23] Feinberg D., ArXiv:cond-mat/0307099 
[24] Prada E. and Sols F., |ArXiv:cond-mat/0307500l 

[25] Lefloch F., Hoffmann C., Sanquer M. and Quirion D., Phys. Rev. Lett, 88 (2002) 197001. 



